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This article gives a brief survey of phase transitions in adsorbed films on 
well defined surfaces of different symmetry of the lattice. The discussion 
concentrates on the effects of periodic changes of the adsorbate — sub- 
strate potential on the structure and thermodynamic properties of ad- 
sorbed films. Different theoretical approaches are briefly reviewed, with 
an emphasis on those which explicitly take into account finite corrugation 
of the surface potential. 


1. INTRODUCTION 


The formation of ordered structures in adsorbed films formed on well de- 
fined crystal surfaces has been attracting a great deal of interest for many years 
[1-5] and the activity in this field is still growing. At present, the literature 
which presents the results of both experimental and theoretical studies is over- 
whelmingly abundant and the achievements are outstanding. The progress in 
this field has gained a particular impetus due to introduction of powerful ex- 
perimental techniques (e.g., neutron scattering [6,7], various forms of electron 
spectroscopy [8-14], light particle, e.g. helium, scattering [15,16], and many 
others [17-21]), the extensive use of computer simulation methods (Monte 
Carlo and molecular dynamics) [5,22-28] and, last but surely not least, owing to 
very intensive theoretical studies [1-3,29-36]. 

The systematic investigation of phase transitions in adsorbed layers began in 
the late sixties [1,37-39] just after the development of the technology for the 
production of the graphite substrate with highly uniform surface [40]. Thomy 
and Duval [37-39] measured adsorption isotherms of simple gases (krypton, 
xenon, methane) on graphite at low temperatures and demonstrated that mono- 
layer films exhibit the two-dimensional counterparts of gas, liquid and solid 
phases. Over the next years, hundreds of papers reporting the new findings were 
published and particularly extensive studies were carried out for the films 
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formed on graphite [41-44]. It soon became clear that several phenomena ob- 
served in adsorbed films do not have simple counterparts in the three- 
dimensional bulk matter. The misfit between the lattice spacing of the substrate 
surface and the preferred lattice spacing of the adsorbate (determined by the 
interactions between the adsorbed particles) leads to the formation of incom- 
mensurate phases, axially ordered phases and other forms of ordering specific to 
surface systems. All those findings clearly demonstrated that spatial nonunifor- 
mities of the gas — solid potential, resulting directly from the lattice structure of 
the crystal surface, are of great importance and often determine the structure of 
adsorbed layers. 

In the case of hexagonal lattices (e.g. graphite surface), the intrinsic lattice 
constant of the adsorbed film, however usually is different from the distance 
between the adjacent minima of the potential wells at different surface cells, but 
nevertheless, the symmetry of the substrate potential and the symmetry of the 
adsorbate preferential structure are the same. Thus, for small misfit, even a 
slight corrugation of the substrate potential may considerably stabilize the 
commensurate phase [1,43], while for large misfit, incommensurate phases are 
formed [45,46]. 

In the cases of square and rectangular lattices, the surface corrugation poten- 
tials exhibit usually much higher periodic variations [1] and there is a natural 
difference in the symmetries of such lattices and the preferred symmetry of the 
adlayers. Thus, the competition between the adsorbate — adsorbate and the ad- 
sorbate — adsorbent interactions must be much stronger and may lead to the 
formation of the new types of ordering in the film. Indeed, several experimental 
studies have revealed the existence of the higher order commensurate phases 
[47,48] and axially ordered phases [12,14,49,50], in addition to the registered 
and incommensurate (floating solid) phases. The mechanisms of phase transi- 
tions between all those different types of surface phases are quite complex and 
are not so well described and understood as in the case of hexagonal surfaces. 

In this paper, we tackle several questions concerning general aspects of the 
formation and properties of the adsorbed films deposited on crystal surfaces of 
different symmetry and our discussion is primarily directed towards the pres- 
entation of recent studies performed for such systems with help of computer 
simulation methods. 


2. THE GAS - SOLID INTERACTION 


Here, we assume that the surface exposed to the gas phase is a perfect single 
plane of a crystal. In such a simple situation we can fully characterize the sur- 
face lattice by its unit lattice vectors & and a, (Figure 1). The interaction of 
the gas atom with the crystal can be represented by the potential v=(T,z), 
where 7 =(x,y) is the two-dimensional vector specifying the location of the 
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adatom in the plane parallel to the surface and z is the distance from the sur- 
face, located at z=0 and assumed to run through the centers of the uppermost 
layer of the crystal atoms. The periodic structure of the crystal surface causes 
that the potential v(T,z) is also a periodic function, such that 

v(£,z) 2 v( - k,á, - k,á,,z) (1) 
where k, and k, are integers. This periodicity property of the potential v(7,z) 
has prompted Steele [51] to represent v(7,z) in the form of the Fourier series 


v, z) = v9(z) + Y v, (z)expliq-7] (2) 


where v,(z) is the interaction potential averaged over the entire surface and the 
sum runs over the nonzero two-dimensional reciprocal lattice vectors q 

q-njb, +n,b, (3) 
where b, and b, are the basic reciprocal lattice vectors and n, and n, are inte- 


gers. 
(100) fee (110) fcc (111) fec 


EA srele E 11 

erel | | Jal | | Te 

EBLJBELGJIENI!J 
E pao] 


3 


8 
i 

3 | 

E i 
C - DP 
ce 

L| 
L| 
d 
>| 
n 


E 

E 

N 

L| 

L| 

|_| a 
BSCE mm 
pela Di] Ble leh RES 


Figure 1. (a) Top view of the fcc(100), (110) and (111) surfaces showing the surface 
unit cell vectors ga, and a,. The lengths of the unit cell vectors for (110) and (111) 


surfaces are expressed in units of a, for the (100) surface. (b) Examples of different 


ordered superstructures on a square lattice. The adsorbate layer unit lattice cell vectors 
are shown for the structure 9 
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The usefulness of the expansion (2) depends primarily on how fast it con- 
verges. It turns out that the convergence of the series is quite fast for the graph- 
ite basal plane and for (100) and (111) faces of the fcc crystal [51,52]. On the 
other hand, for the (110) plane of the fcc crystal it is necessary to include sig- 
nificantly more terms in order to obtain reliable results [53]. 

The expressions for the Fourier coefficients derived by Steele have rather 
limited application to real systems, since the Lennard-Jones potential is not the 
best choice for the representation of two body interactions in adsorption on 
ionic crystals, e.g., NaCl and MgO, or on metals. 

Recent experimental [54] and theoretical [55,56] studies have demonstrated 
that even the closed packed faces of metal crystals, e.g. Pt(111), show quite 
strong periodicity of the surface potential, against a wide-spread belief [57,58] 
that such surfaces can be considered as essentially noncorrugated. ` 

In fact, there are not so many direct and reliable estimations of the surface 
corrugation potential, apart from rather well known examples of simple gases 
on graphite [59], lamellar dihalides [60], and on some close packed metal sur- 
faces [16]. It should be noted that even in the case of simple gases adsorbed on 
graphite, for which the interaction potentials are believed to be known quite 
well, it turns out that the corrugation effects are probably higher than obtained 
from theoretical calculations. It has been demonstrated by computer simulations 
[61-66] that in order to obtain a good agreement with experimental data the 
periodic variations of the surface potential need to be higher. To take this effect 
into account, the expression (2) has been modified [61] by introducing the ad- 
justable "corrugation parameter" V, , so that eqn. (2) takes the following form: 


v(E, z) 2w(z)4- v, v, (zexpliq .£] (4) 


For the values of V, >1, the effects due to the corrugation are enhanced with 
respect to those predicted by eqn. (2), while for V, «1 they are weakened. 

In many cases, to study basic effects of the surface potential corrugation and 
the surface lattice symmetry it suffices to retain only very few leading terms in 
the expansion (2) (or (4)) and consider the Fourier amplitudes as adjustable 
parameters. A good example of such approach is the study performed by Bruch 
and Venables [67], who considered the conditions that must be satisfied by the 
system geometry and the binding potential in order to stabilise uniaxially or- 
dered structures in monolayer films. The above assumption is justified by the 
known property of the Fourier coefficients v.(z), which usually decay rapidly 
with || and z. In many cases this decay is exponential [51], i.e., 


vy)  expl-|q|z] (5) 
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From the above property of the corrugation potential, we immediately con- 
clude that its direct effect on the admolecules located beyond the first layer is 
bound to be very weak, if any at all. It does not mean, of course, that the struc- 
ture of the second and higher layers is unaffected by the properties of the sur- 
face potential. The attraction between admolecules and the surface is responsi- 
ble for the formation of subsequent adsorbed layers and its strength is a chief 
factor determining the film growth mode [34], as well as internal structure of 
particular layers. The effects due to the surface corrugation are transmitted to 
the higher layers through the interaction with the admolecules from the first 
layer, which are under a strong influence of the surface corrugation potential. 
These effects have been demonstrated by the molecular dynamics simulation of 
Phillips and Shrimpton [66], as well as by the Monte Carlo study of Patrykiejew 
et al. [68,69]. 

The properties and structure of adsorbed films are controlled by the com- 
bined effects due to admolecule — admolecule and admolecule — substrate inter- 
actions as well as by the imposed thermodynamic constraints (specified tem- 
perature, film density, bulk gas pressure). Of course, in real situations, as well 
as in computer simulations, the effects due to finite size of the surface may con- 
siderably affect the structure of adlayers [68], but we shall not discuss such 
problems here. 

In general, at low temperatures the adsorbate — adsorbate interaction tends to 
enforce the formation of hexagonal close-packed (hcp) solid phase in mono- 
layer films, while the corrugated admolecule — substrate potential favours the 
formation of registered structures. A simple criterion enabling to classify the 
adsorbed films with respect to their structure has been developed by Park and 
Madden [70]. Assuming that the substrate surface lattice unit cell vectors are a, 
and d, and the overlayer unit cell is characterized by the unit vectors é, and e, 
the relation between these two sets can be written as 


€; a, Ay jJ a, 
Park and Madden have classified the adsorbed layers with respect to the be- 
havior of detlaz, | and singled out three cases. The first case groups all systems 


for which detlar, | is an integer number. This situation corresponds to the regis- 


tered or commensurate adsorbed films. The second case involves the situations 
in which det]o;, | is a rational number. When this condition is satisfied the ad- 
sorbed film forms the so-called “high-order” commensurate phase, with only a 
certain fraction of adatoms located directly over the adsorption sites. Finally, 
when detlo;, | is an irrational number (case three) the adsorbed layer is incom- 
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mensurate with the substrate surface lattice. The above classification parame- 
trises the adsorbed layer structure with respect to the relative sizes of the surface 
and adlayer unit cells. The area of the surface unit lattice cell is given by 
à, Xd, , and the area of the adsorbed layer unit cell is equal to e, xe, , so that 


detlo;, | is equal to the ratio of these two areas. 


The criterion of Park and Madden vastly oversimplifies the reality. The dis- 
tinction between the incommensurate and the high-order commensurate phases 
may cause problems, as it is usually possible to approximate irrational number 
by a suitably chosen rational number that falls into the region of experimental 
resolution. Aubry [71] has even proposed to consider any incommensurate ad- 
sorbed layer as sufficiently high-order commensurate phase. When the condi- 
tions change and the density of the film increases or decreases, it is supposed 
that the resulting changes in the film structure are due to a series of first-order 
transitions between different high-order commensurate phases. This leads to the 
so-called "devil's staircase" of phase transitions. 

Another weakness of the above criterion is that it does not include incom- 
mensurate phases exhibiting domain walls networks [45,72]. In such cases the 
incommensurate phase is composed of large commensurate domains separated 
by walls which can have different structure, thickness and orientation. 

A useful measure of the relative size of the adsorbate and the surface lattice 
unit cell is the so called dimensional incompatibility parameter [60], defined as 


I 2 (a-n)!n (7) 


where a is the surface lattice constant and r, is the distance between adjacent 


adsorbate atoms in the surface phase. Of course, in the case of a rectangular 
symmetry of the surface lattice, one can define two different dimensional in- 


compatibility parameters /, and/,, by taking a=|a,| and a=|ā,|. The behav- 
ior of films formed on corrugated surfaces is closely related to the magnitude of 
the parameter 7 [60,73]. 


3. LATTICE GAS MODELS 


The interaction between the particles adsorbed on a lattice can be repre- 
sented by the potential u(r), where r assumes discrete values: 


rz4Kk'a? lai (8) 
In the above a, and a, are the lengths of the unit cell vectors and k and / are 


integers. The possible ordered structures that can appear in the system are pri- 
marily determined by the lattice symmetry and the properties of the potential 
u(r). Usually, u(r) is a sum of contributions arising from repulsive and attrac- 
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tive interactions. For example, assuming that u(r) is given by the Lennard- 
Jones potential, the effect of repulsive forces can be related to the relative size 
of adatoms, measured by the parameter o and the size of the surface unit cell, 
measured by the lattice constants a, and a,. It is convenient to introduce the 
reduced quantities and express all the distances in units of, say, a, . Note that in 
the case of a square lattice, eg. (100) plane of the fcc crystal, a, = a, , while for 


the (110) plane of the fcc crystal a, = «2a, . Also, the energy-like quantities, as 
well as the temperature, can be expressed in units of €. 

Lattice gas model provides a simple explanation for the experimentally ob- 
served changes in the critical temperature of adsorbed monolayers resulting 
from the changes of the dimensional incompatibility between adsorbate and 
adsorbent (cf. eqn. (7)). Note that in the case of a lattice gas model, a simple 
mean field approximation gives 


U 
TUI m sum 9 
P TUE (9) 
where 
U „m 7-05 Y u(r,) (10) 
j 


where the sum runs over different shells of neighbors. In the particular case of 
Lennard-Jones interaction potential, truncated at 2.50", U „m is given by [73] 


U,,,, =-8e[C,,0°" - C,0* (11) 


where C,, =1.015997141 and C, =1.156625. Thus, the critical temperature is 
proportional to U,,,,, which varies with o`. Of course, the same is true for 
more rigorous treatments, as it is illustrated by the results shown in Figure 2a 
and 2b. The only available experimental data that clearly demonstrate changes 
of the two-dimensional critical temperature with dimensional incompatibility 
correspond to adsorption on a honeycomb lattice of lamellar dihalides [60] (cf. 
Figure 2c). Nevertheless, the qualitative picture should be the same in the case 
of adsorption on a square lattice. Of course, the predictions stemming from lat- 
tice gas models do not directly apply to adsorption on rather weakly corrugated 
surfaces of lamellar dihalides and are, at the most, only qualitatively correct. 
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Figure 2. Two-dimensional critical temperature versus U,,,, for the Lennard-Jones parti- 
cles of different diameter, o` =o /a, (shown in the figure) on a square lattice, obtained 
from mean field approximation (part a) and from the coherent anomaly theory (part b). 
Part c shows the ratio of two-dimensional and three-dimensional (bulk) critical tem- 
peratures for simple gases adsorbed on lamellar dihalides (from Ref. 73) 


Lattice gas models work much better for strongly adsorbed (chemisorbed) 
species. Good examples of such a system are the monolayers of atomic hydro- 
gen adsorbed on Pd(100) [11] and on W(100) [74], Cl adsorbed on Ag(100) 
[75] and O atoms adsorbed on Cu(100) [76]. It is particularly interesting that 
even very simple lattice models allow to obtain rather good agreement with 
experimental data. For instance, hydrogen adsorbed on W(100) orders into 
c(2x2) phase at low temperatures. The simplest possible model which predicts 
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that form of ordering is the Ising-like model with nearest-neighbor repulsive 
interactions only. Monte Carlo simulation for such a model [77] led to rather 
good agreement between the calculated and measured temperature changes of 
the intensity of a LEED spot for the c(2x2) phase. 
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Figure 3. Coverage dependence of the (1/2,1/2) LEED beem height at 300K for the 
atomic Cl adsorbed on Ag(100) surface (part a), and the structure factor obtained from 
Monte Carlo simulation on the 72x72 lattice (part b). (Adapted from Ref. 75) 


A quite similar model was also found suitable to describe the dissociative 
adsorption of chlorine on Ag(100) surface. A monolayer film of atomic Cl on 
Ag(100) was found to form the c(2x2) phase [77,78]. Taylor et al. [75] per- 
formed Monte Carlo simulations for a model with infinite nearest-neighbor 
repulsion on a 72x72 lattice and found a good agreement with experimental 
data. In particular, the calculated structure factor versus coverage reproduces 
quite well the measured LEED beem height (see Figure 3). A rather convincing 
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proof of adequacy of the Ising model representation for that system is the 
“Fisher renormalized” exponent 8 /(1—a@) = 0.12 evaluated from experimental 
data and which agrees very well with the exact Ising value of 0.125. 

Although the formation of c(2x2) ordered phase is quite common in atomic 
monolayers on metals, there are also known numerous examples of systems 
which also show other ordered phases. For example, atomic oxygen adsorbed 
on Rh(100) was found to form the p(2x2) structure of the density 0.25 [79], 
which then undergoes transition to a denser c(2X2) structure. The same was 
found for O/Pd(100) [67] and Se/Ni(100) [80]. The phase diagram for the last 
system was evaluated by reflection high-energy electron diffraction (RHEED) 
[80] and was also studied by Monte Carlo simulation [80,81]. Analysis of ex- 
perimental data, model calculations and symmetry arguments led to the conclu- 
sion that Se/Ni(100) phase diagram is topologically equivalent to the Ashkin- 
Teller model [82]. In this model, each lattice site is characterized by two Ising 
spins, s, and 1, , and the Hamiltonian reads: 

Hz-JY (sss, * tt) - AY sts t, (12) 

«i, j> <i,j> 

Critical properties of that model are believed to be non-universal and dependent 
on the ratio A/J . Monte Carlo study was performed [80] for a simple lattice 
gas model with interactions up to fourth nearest neighbor which yields the 
phase diagram belonging to the universality class of Ashkin-Teller model. Fig- 
ure 4 shows that the experimental Se/Ni(100) phase diagram and the phase dia- 
gram obtained for the above lattice gas (Ashkin-Teller-like) model are quite 
similar, indeed. 

Extension of the lattice gas formalism from two- to three-dimensional ad- 
sorption system is straightforward [34,83]. Assuming that the adsorbate — sub- 
strate interaction is represented by the potential v(z), where z is the distance 
from the surface, the Hamiltonian for such a model reads 


1 

H-Y v(z,)n, t2, mna, (13) 
i ij 

Note that in the lattice gas language, z, can take upon only discrete values de- 

termined by the distance of the i-th site from the surface, given by the number 

of layer |. Depending on the properties of the potentials v(z) and u(r) the 

system shows different behavior [34,84-89]. 
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Figure 4. A comparison of experimental and theoretical phase diagrams for selenium on 
Ni(100) surface. (From Ref. 80) 


Systematic analysis of lattice gas model representation of multilayer adsorp- 
tion was done by Pandit et al. [34] for the potential v(z) given by 


v(I) =V,6,, + BI? (14) 


where V, is the energy of adsorption for the first adsorbed layer. They pre- 
sented a very elegant and detailed discussion of possible scenarios for the film 
growth on the basis of mean field theory for lattice gas model. More recently 
this problem was also considered by Prasad and Weichman [90], who used both 


the mean field theory and the renormalization group arguments. 
Here we briefly discuss the main predictions stemming from the lattice gas 


model assuming that B =V, in eqn. (14) and neglecting all but the first nearest 
neighbor interactions between adsorbate atoms. 
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Figure 5. An example of experimental adsorption isotherm for Argon on MgO at 


T = 60.8 K (adapted from Ref. 91) (part a) and isotherms obtained from Monte Carlo 
study (part b) for a cubic lattice gas model (from Ref. 88) 


When the ratio V,/u (u=u(ryy)) is sufficiently high, i.e., when the surface 


is highly attractive towards the adsorbate, the film grows in a layer-by-layer 
mode (see Figure 5) which was observed in many real adsorption systems 
[3,41,60,91]. At low temperatures it corresponds to the presence of the series of 
layering transitions which terminate at the corresponding critical points T, (l). 


From the mean-field theory it follows that 

limT, (1) = T7" (15) 
where T?” is the critical temperature of the uniform bulk gas, while the exact 
result is [34] 

limT, (I) 2T, (16) 
where T, is the roughening temperature, which in general is much lower than 


T?” and only slightly exceeds T, (1) and the following relation is predicted [92] 
to be obeyed 


T, - T. (I) œ (Inl)*, | (17) 


In the above situation the adsorbate wets the surface at any temperature down to 
T 20. 
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As the adsorbate — substrate interaction becomes weaker and the following 
condition is met 


= -l 
Vol ba (18) 
u Isl 

the adsorbate does not wet the surface at T =0 , but may still exhibit complete 
wetting at the temperatures which exceed the wetting temperature 7,,. Thus, 
below T, the film thickness remains finite up to the bulk condensation point, 
while for T >T, we observe again the formation of macroscopically thick ad- 
sorbed film as the chemical potential approaches its bulk coexistence value. 
That change of the film behavior at T, marks the so-called wetting transition 
which can be either the first-order or the continuous transition (critical wetting) 
depending on the value of V,/u and on the range of molecular interactions in 
the system [93]. When the wetting transition is first-order it is often preceded by 
the prewetting transition between the thin and thick film. The prewetting tran- 
sition occurs at the chemical potential below the bulk coexistence value. As the 
temperature increases the prewetting transition terminates at the corresponding 
critical point. Due to the finite extent of thin and thick films, on the both sides 
of the prewetting transition, it belongs to the universality class of the two- 
dimensional Ising model. This situation corresponds to the so-called intermedi- 
ate substrate regime. One can interpret the wetting transition at 7, as an inter- 
face unbinding transition [94]. At temperatures below T,, in the non-wet re- 
gion, the gas-liquid interface is very tightly bound to the substrate surface and 
only a small density enhancement near the wall occurs. In the wet state this 
density enhancement is large and even exceeds the density of the bulk liquid 
over quite large distances from the surface and the interface between the liquid 
layer and the bulk gas phase is a smooth, delocalized object of large height am- 
plitudes z,(x, y). The wetting transition temperature grows as the ratio V, /u 
becomes lower. When it reaches a certain limiting value, the wetting tempera- 
ture reaches the bulk critical point of the adsorbate, and for still lower V, /u one 
enters the so-called weak substrate regime. In this case the adsorbate does not 
wet the substrate under any conditions and adsorption is always finite. 
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Figure 6. (a) Phase diagram in the coverage (0) — temperature (T° =kT/J where J is 
the exchange constant between nearest neighbors) plane of a nearest neighbor lattice gas 
model an a simple cubic lattice with a free surface, and a potential (14) with 
B - V, 2 2.5J . (b) The corresponding phase diagram in the grand-canonical ensemble 
u-(u-p,)/J where y, is the chemical potential at the bulk gas — liquid coexistence) 
Note that the adsorbate — substrate potential was cut off above the fourth layer, so that 
the curve for the layering transitions in the fifth and higher layers merge at the bulk 
coexistence curve (from Ref. 88) 


Within the strong substrate regime the properties of adsorption systems ap- 
pear to be very sensitive to even small changes in the relative strength of adsor- 
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bate — adsorbate and adsorbate — substrate interactions. Figures 6 and 7 show 
examples of phase diagrams resulting from Monte Carlo study of a lattice gas 
model with nearest neighbor adsorbate — adsorbate interaction but long range 
surface potential (given by eqn. (14) with B —V,) [88]. In the case of strong 
substrate potential the film grows via a series of simple layering transitions in- 
volving condensation in one layer each. For a weaker surface potential one 
finds, however, a different sequence of layering transitions in which the first 
three layers condense together at low temperatures. This system shows also an 
example of a surface triple point at which the layering transitions of layers 1+2 
and layer 3 coexist. Note that the critical point of the layering transition for 
layers 1+2 is considerably higher than the corresponding critical points for 
higher layers. This results from the enhanced effects of adsorbate — adsorbate 
interaction in the condensed two surface layers as compared with a single layer. 
Of course, even for a transition involving simultaneous condensation in several, 
but finite, number of layers the critical behavior corresponds to the universality 
class of a two-dimensional model, as the correlation length in the direction per- 
pendicular to the surface stays finite. 
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Figure 7. The same as in Figure 4.15 but for a weaker substrate potential with 
V, =0.93J (from Ref. 88) 
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A sufficient decrease of substrate potential finally leads to the crossover 
from the strong to intermediate substrate regime. This crossover may be also 
induced by adding further neighbor (attractive) interactions between adsorbate 
particles [89]. Figure 8 shows a comparison of the phase diagram obtained for a 
simple model with the first nearest neighbor interaction with the phase diagrams 
corresponding to the model with added interaction between the third nearest 
neighbors. 


4. FILMS ON SURFACES WITH FINITE CORRUGATION 


As it was already demonstrated in section 2, atoms adsorbed on crystalline 
surfaces experience corrugated surface potentials. In the case of physisorption 
periodic variations of the gas — solid potential are usually small as compared 
with the adsorption energy and often are also smaller than or comparable with 
the thermal energy of the adsorbed atoms. In such situations the lattice gas 
models discussed in the previous section do not describe well the properties of 
adsorbed layers. Off lattice movement as well as out-of-plane excitations of 
adsorbed atoms give rise to many new phenomena that can not be included into 
the lattice gas formalism. 


Figure 8. A comparison of the phase diagram for the lattice gas models with the first 
nearest neighbor interaction (part a) with the phase diagrams for the model with added 
second neighbor attractive interactions (parts b and c). Parts b and c show the results for 
R=u,,,/u,,=0.05 and 0.08, respectively. Phase diagram in part a was taken from Ref. 


88, while those shown in parts b and c from Ref. 89 
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In the case when the gas — solid potential exhibits periodic variations due to 
the lattice structure of the surface, a competition between the surface corruga- 
tion and the adsorbate — adsorbate interactions becomes a major factor deter- 
mining the structure of adsorbed layers. In general, the adsorbate — adsorbate 
interaction tends to enforce the formation of the hexagonal close-packed (hcp) 
incommensurate solid phase in the monolayer film at low temperatures. On the 
other hand, the corrugated surface potential favors the formation of registered 
(or commensurate) structures. The actual structure of the film may be quite 
complicated and characterized by highly nonuniform spatial distribution of ad- 
sorbed atoms within the film. In the case of the most intensively studied ad- 
sorption on substrates with a hexagonal symmetry of the surface lattice, such a 
graphite [41,95-97], boron nitride [98,99], dense (111) plane of metal crystals 
(Pt, Cu, Ag) [13,100-103] and lamellar dihalides [60], both the commensurate 
and the incommensurate solid phases have the same symmetry and differ only 
by the spacing between neighboring atoms and by the orientation of the ad- 
sorbed layer lattice relative to the surface lattice [13,31,104-106]. In many cases 
the structure of the incommensurate phase can be described by the domain wall 
formalism [45,107]. 

In the case of adsorption on surfaces of a square and rectangular symmetry, 
such as low index planes {(100) and (110)} of the face centered crystals of 
various metals [14,108-112] or on the (100) plane of ionic regular crystals of 
MgO [91,113-117] and alkali halides [118-120], the situation looks quite differ- 
ent. Competing adsorbate — adsorbate and adsorbate — substrate interactions 
give rise to the formation of various registered superstructures [10], uniaxially 
ordered phases [121] as well as to various incommensurate phases [122]. 

Several important informations about the ordering in two-dimensional films 
formed on crystalline surfaces can be derived from the ground state calculations 
[52,67]. Bruch and Venables [67] considered the relation between the surface 
lattice and the adsorbed layer geometry in two-dimensional films adsorbed on 
crystals of different symmetry. They concentrated the discussion on the effects 
due to the corrugation potential and formulated explicit conditions for uniaxial 
registry. In that treatment the adsorbed film was assumed to be uniform in a 
sense that it was strictly two-dimensional and possessed a well defined lattice 
structure. The results obtained by Bruch and Venables must be considered as a 
sort of zeroth-order approximation only, since they did not attempt to minimize 
the system total energy, but rather focused on the calculations of the registry 
energy exerted by the surface corrugation potential only. 

Some simple examples of the ground state calculations were reported by Pa- 
trykiejew et al. [52,53,69] in the case of adsorption on model (100) and (110) 
planes of atomic face centered cubic crystals, assuming that the adsorbed film 
forms a regular hexagonal lattice or simple registered phases such as (1x1), 


c(2x2) and (2x1). It was assumed that the adsorbate — adsorbate interaction 


266 A. Patrykiejew 


potential u(r) is given by the Lennard-Jones function while v,,(7) was repre- 
sented as 


Vap (4) = min[v(z, z)] (19) 


with v(T,z) given by eq. (4). 

In the case of an incommensurate hexagonal phase consisting of N atoms, 
the energy corresponding to a given distance between the first nearest neighbors 
(œ) and to an also specified orientation of the adsorbed film with respect to the 
surface lattice (epitaxial rotation) [104-106] represented by the angle 2 one has 


e'(0;,9)- 1alce (c^ /a;) -ct(o* a Y] sen | (20) 


where all energies are expressed in units of 2,, and all lengths in units of œ. 


This energy must be then minimized with respect to the both oz and & in order 


to find the stable configuration. One readily observes that in the thermodynamic 
limit (N — e) such a perfect, incommensurate, hexagonal phase has no pre- 


ferred orientation [106] and the energy is given by 


e' (o1) -1ijc*(s 7o; }? LE (o la; |+ vipo (21) 


E 


where v, , is the zeroth-order Fourier coefficient of v; (c ) given by 
` 1 TEEN 
Vap,o -—[ V2p G Az (22) 
a, ? 


where the integration is performed over the surface lattice cell of area a,. The 
parameters C; and C; can be readily evaluated and their values depend on the 


assumed range of molecular interaction [123]. 
In the case of registered phases, in which all adsorbed atoms are located di- 


rectly above the surface potential minima (7; = (0.5,0.5a;)} we have 
e (7) -8c;o"" - c;o* ]+ vi (@7) Q3) 


where the magnitudes of the parameters Cj, and C; depend on the assumed 
commensurate structure r [69]. 

Using such an approach the stability regions corresponding to the hcp and 
the registered (1x1), c(2x2) (and (2x1) in the case of adsorption on the (110) 
plane of an fcc crystal) phases in adsorbed films formed on surfaces character- 
ized by different size of the adsorbate atoms, given by o”, and different corru- 
gation of the surface potential have been estimated (e.q., see Figure 9). 
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Of course, such calculations do not provide any exact solution to the problem of 
ground state properties of incommensurate phases as they ignore lattice distor- 
tions which arise from defects and local elastic relaxations, as discussed by 
Vives and Lindgard [106] as well as from the formation of domain-wall net- 
work structure [45]. 

To prepare a ground for a full appreciation of the surface corrugation effects 
on the finite temperature behavior of films adsorbed on crystalline surfaces, we 
now briefly discuss the properties of films formed on a “flat” noncorrugated 
surface. When the surface potential possesses a deep minimum at a certain pre- 


ferred distance from the surface, z,, then the adsorbed film is essentially two- 
dimensional at low temperatures. Deviations from planarity of adsorbed layer 


may be due to temperature rise [124,125] as well as may result from the forma- 
tion of the second and higher layers [68]. 


8 9 1.0 1.1 1.2 


c 
Figure 9. Regions of stability of different surface phases for the Lennard-Jones particles 
adsorbed on the (100) plane of model fcc crystals deduced from the ground state calcu- 
lations for the systems of the number density equal 1.0 (from Ref. 52) 


The behavior of strictly two-dimensional Lennard-Jones systems has been 
intensively studied by various theoretical approaches [126-129] as well as by 
computer simulations [25,127,130-134]. In general, the phase diagram is quite 
well known. It has been found that it exhibits two-dimensional counterparts of 
all familiar states of matter in three-dimensional space: gas, liquid and solid (see 
Figure 9). The triple point temperature has been estimated quite accurately as 


equal to 7; - kT/€ 20.40 € 0.015 [127,134]. Existing estimations of the criti- 
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cal temperature for that system are much less conclusive [127,135] and obtained 
values range from about 0.5 to 0.56. 

Another controversy about the properties of strictly two-dimensional Len- 
nard-Jones system concerns the mechanism of melting. Already in the 1930's 
Peierls [136] argued that truly long-ranged positional order can not exist in two- 
dimensional systems, so that two-dimensional crystalline order is impossible. 
This has been proved later by Mermin and Wagner [137]. In 1973 Kosterlitz 
and Thouless [29] proposed a theory dislocation-mediated melting for a two- 
dimensional systems. That theory has been later developed further by Halperin 
and Nelson [30,138] and by Young [139]. From those theoretical works it fol- 
lows that two-dimensional systems possess only quasi-long-renged positional 
order characterized by algebraic decay of the two-particle correlation function. 
The KTHNY theory predicts that melting in two-dimensions occurs via two 
continuous phase transitions. In the first transition, due to dissociation of dislo- 
cation pairs, the system looses the quasi-long-range positional order, but retains 
the quasi-long-range orientational order. This corresponds to the formation of 
the so-called hexatic phase. The second stage of melting is connected with the 
disclination-unbinding transition, occuring at higher temperature and leading to 
the formation of isotropic two-dimensional liquid. 
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Figure 10. Schematic representation of the phase diagram for the two-dimensional Len- 
nard-Jones system 


The scenario predicted by the KTHNY theory is considerably altered by the 
presence of periodic substrate potential and the actual mechanism of melting 
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depends on the symmetry and the size of the substrate surface as well as on the 
amplitude of the periodic surface potential [138]. In general, surface corrugation 
is expected to wash-out the disclination-unbinding transition, and hence the 
melting becomes a one-step process. This transition may be continuous, as for 
the flat surface, or first-order. The only exception is the melting of a hexagonal 
adsorbate film on a square lattice characterized by weak corrugation of the ad- 
sorbate — substrate potential. In this case the theory proposed by Nelson and 
Halperin [30,138] predicts that the disclination-unbinding transition is replaced 
by an Ising-like transition. Thus, one expects to observe two liquid-like phases 
with different orientational symmetry. 

There are also other theories [140-143] that predict the usual first-order 
melting in two-dimensional systems. The concept of first-order melting was 
also advocated by Abraham [25,144], who argued that the loss of positional 
order is negligible in real systems due to their limited, though macroscopic, 
size. 
Computer simulation studies [25,133,134,144-146] as well as experimental 
data [147-151] do not provide an univocal answer to the problem of melting in 
two-dimensions. In real adsorption experiments the assumption of two- 
dimensionality of adsorbed film is seldom fulfilled. Out-of-plane moyement of 
adsorbed particles is expected to be important factor which influences the prop- 
erties and stability of all possible phases in monolayer films [68,123]. 

Monte Carlo simulations performed for a series of two-dimensional systems 
of Lennard-Jones particles with o° = 0.9 on surfaces with different corrugation 
[152] demonstrated that the mechanism of melting on corrugared surfaces is 
indeed considerably different from that predicted for a flat surface. 

In order to determine the order and nature of the melting transition the be- 
havior of the suitable bond-orientational order parameters [145] and their 
fourth-order cumulants for systems of different size were studied. The adlayers 
with a four-fold symmetry, as in the case of the both (1x1) and c(2x2) phases 


formed on the (100) plane of an fcc crystal, and the layers with a six-fold sym- 
metry can be characterized by the following bond-orientational order parame- 
ters [152,153] 


v, =H Ye] (24) 
N, lo otl 


In order to detect axially ordered hexagonal phase one can define the bond- 
orientational order parameter of the form 


y= s TE solo (25) 
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In the above the first sum runs over the all molecules in the system, while the 
second sum is taken over the first shell of neighbors, N, is the total number of 
"bonds" in the system and @ is the angle between the bond and the reference 
axis, taken to be the x-axis of the surface lattice. 
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Figure 11. Bond-orientational order parameters y, (part a) and y, (part b). The results 


shown in part (a) correspond to the systems with c^ =1.20, three different values of the 
corrugation parameter y, and three different system sizes (shown in the figure) while 


the results shown in part (b) were obtained for y, 20.6 and three different system sizes 
(also shown in the figure) (from Ref. 153) 


When the adsorbed layer is in a disordered fluid state, all the above defined 
bond-orientational order parameters should be equal to zero since all possible 
mutual orientations of “bonds” appear in the system with the same probability. 
In the case of registered phases of a square symmetry, y, is expected to be 
equal to unity, while the remaining order parameters y, and y; should be 
equal to zero. In the hexagonal close packed phase y, should approach zero 
while y, should be equal to unity. The behavior of y; depends on the orienta- 
tion of the adsorbed layer with respect to the surface lattice and may assume 
different values. When the film exhibits uniaxial registry but retains perfect 
hexagonal structure than y; —1.0. On the other hand, when the hexagonal ar- 
rangement is not perfect, we expect nontrivial behavior of y; , as well as of y, 
and Y, . Examples of the changes of y, and y, with temperature for adsorbed 
films formed on highly corrugated systems of different size are shown in Figure 
11. It is quite clear that in all cases the adsorbed layers exhibit square symmetry 
at low temperatures and undergo rather abrupt, though continuous, disordering 
as the temperature increases. The absence of any finite size effects is a clear 
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evidence that no longwavelength fluctuations develop. The finite size effects are 
only manifested in the behavior of residual y, (cf. Figure 11b), but this merely 
reflects gradual approach towards the thermodynamic limit. 

It is well known [154,155] that in the case of a second order phase transition 
the fourth-order cumulants of the order parameter reach a trivial limit of 2/3 at 
the temperatures well below the transition point and zero at the temperatures 
well above the transition point. At the transition temperature the cumulants for 


different system sizes should reach a common nontrivial intersection point U" 


with the magnitude of U” depending on the universality class of the phase tran- 
sition. The results presented in figure 12 seem to confirm the predictions of the 


KTHNY theory. At the temperature T’ ~0.4 the bond-orientational order pa- 
rameters y, and y; both exhibit a sudden, though small, drop. 
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Figure 12. Bond-orientational order parameters y, and y; for the two-dimensional 
system of 894 particles with o* =0.9 adsorbed on the surface of low corrugation, char- 
acterized by V, =0.1, and the size of 32x28 surface unit cells as well as the fourth-order 
cumulants U, (L) for different system sizes obtained from Monte Carlo simulation (from 
Ref. 152) 


Thus, above that temperature the system retains a considerable orientational 
order. According to the KTHNY theory this first transition, due to dissociation 
of dislocation pairs, transforms the solid phase into the hexatic phase. Then a 
second transition occurs at the temperature of about 0.485, which leads to the 
loss of orientational order in the system. The behavior of U, , also depicted in 


Figure 12, shows that the cumulants for different systems sizes do show a 
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common intersection at U” close to 0.61. This value of the cumulant fixed point 
correspond to the universality class of the Ising model [154]. 

Real adsorption systems usually show some effects of out-of-plane motion 
of adsorbed particles, even at low temperatures. The increase of the pressure, or 
the chemical potential, for the monolayer film may induce many new phenom- 
ena. What actually happens to the film depends crucially on the strength of the 
adsorbate — adsorbate and the adsorbate — adsorbent interaction, as already dis- 
cussed in sec. 2. Also the surface corrugation is expected to be a very important 
factor influencing the film behavior. 

Possible scenarios may include the appearance of layering transitions, for- 
mation of compact three-dimensional crystallites, promotion of the second layer 
due to melting and disordering of the first adsorbed layer as well as simple de- 
sorption of particles due to increase in their kinetic energy. Also, the structure 
of adsorbed layers adjacent to the substrate surface may change upon the for- 
mation of higher layers, so that various commensurate — incommensurate tran- 
sitions may occur in the film. 

The problem of monolayer stability is closely related to the wetting of a sub- 
strate by an adsorbed film. In the resent paper by Phillips [156]. is was demon- 
strated via the molecular dynamics simulation and analytic solution of the elas- 
tic continuum models that three different mechanisms may be responsible for an 
instability of the adsorbed monolayer and the promotion of the second, as well 
as higher, layers. In a strongly adsorbed film, the promotion of the second layer, 
under the specified thermodynamic conditions of temperature and pressure, 
occurs mostly due to migration of the atoms located at the edges of adsorbed 
islands of a compact two-dimensional solid layer. This process may be also 
accompanied by, much more rarely occuring, movement of individual atoms 
from the interior of the island to the second layer, if they happen to gain enough 
kinetic energy. These two mechanisms were observed also in the case of melt- 
ing in dense monolayer films as will be discussed later. The third possible 
mechanism leading to the promotion of higher layers may be a sudden promo- 
tion of groups of adsorbed atoms due to the response to structural instability in 
the highly compressed solid layer. This last mechanism is most likely to occur 
in dense and relatively weakly adsorbed films in which stresses and strains are 
of particular importance. 

Of course, the above mentioned effects are most likely to occur in adsorbed 
films formed on rather weakly corrugated surfaces, since such systems are most 
likely to form islands of a high density solid phase. In the case of strongly cor- 
rugated surfaces, on which the registered films are usually found, the promotion 
of the higher layers is not likely to occur. It can be rather expected that only a 
sufficient increase in the temperature may induce disordering which is accom- 
panied and/or followed by desorption. 
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The stability of the registered structures in the first layer appears to be quite 
sensitive to even small changes of the surface corrugation as well as to the den- 
sity changes in the film, as it was recently demonstrated by Patrykiejew et al. 
[69]. An interesting phenomenon observed in adsorbed films studied by Monte 
Carlo simulation is the reconstruction of the first adsorbed layer due to the for- 
mation of the second layer. Namely, it was found that, in some cases, the dense 
monolayer films formed on a sguare lattice assume incommensurate structure of 
hexagonal symmetry, while the appearance of the second layer causes a gradual 
restructuring of the first layer, which recovers the ordered epitaxial structure of 
a square symmetry (see Figure 13). Thus, the appearance of the second layer 
exerts a force which pushed the adsorbed atoms from the first layer back into 
the registry positions. 


Figure 13. Snap shot configurations for the adsorbed layer of Lennard-Jones particles of 
o` =0.8 on the (100) plane of the model fcc crystal with the corrugation parameter 
V, =0.9 recorded at the temperature 7° =0.02 and the total film number density equal to 
p, = 1.5625 (part a) and 2.0 (part b). Particles located in the first layer are represented by 
open circles while those from the second layer by filled circles (from Ref. 69) 
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Figure 14 shows the phase diagram derived from the Monte Carlo simulation 
for another system, which exhibits the registered (1x1) and the incommensu- 


rate structures in the monolayer region and also incommensurate second layer. 
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Figure 14. Phase diagram for the Lennard-Jones adsorbed bilayer film formed on the 
(100) plane of fcc crystal characterized by the corrugation parameter y, =0.8. Part a 
shows the phase diagram in the plane (p,,7°). Filled points are the results of grand 
canonical Monte Carlo calculations, while stars and diamonds correspond to the canoni- 
cal ensemble calculations. The triple line of the solid — liquid — gas coexistence in the 
second layer is shown as a broken vertical line. Part b shows the phase diagram in the 
(uU, T^) plane derived from the grand canonical ensemble Monte Carlo simulation. The 
gas — commensurate (1x1) phase transition points are shown as filled points, the com- 
mensurate — incommensurate transition points in the first layer are represented by open 
points, while the layering transition points in the second layer are marked by open dia- 
monds. The location of the gas — (1x1) phase transition at the zero temperature was 
obtained from the ground state calculations, while the locations of the commensurate — 


incommensurate and layering transitions at the zero temperature were estimated by 
extrapolation of the results obtained for finite temperatures. (From Ref. 68) 
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Figure 15. The heat capacity curves for the Lennard-Jones adsorbed films of different 
density (shown in the figure) formed on the (100) plane of fcc crystal characterized by 
the corrugation parameter V, =0.8. The peaks at the temperature of about T° =0.4 ob- 
tained for the systems with the two highest values of the density correspond to the tri- 
ple-point melting of the second layer, while the remaining peaks result from the melting 
of the monolayer incommensurate solid phase (from Ref. 68) 


The data shown in Figure 15 demonstrate that the temperature of the melting 
transition of the incommensurate monolayer solid phase moves towards higher 
values when the density of that phase increases. Only for the films which ex- 
hibit the formation of the second layer the melting temperature is constant and 
equal to about 0.5. Quite similar effect was found experimentally for nitrogen 
adsorbed on graphite [44]. It was demonstrated that it is connected with a dif- 
ferent mechanism of the melting transition in submonolayer and dense mono- 
layer solid films. Submonolayer solid does not occupy the entire surface but 
rather forms high density islands which coexist with a dilute phase. Upon melt- 
ing the more or less uniform fluid phase is formed and spreads over the entire 
surface. This can not happen in already dense monolayer since there is not 
enough space for the fluid phase. The only possibility for the liquid monolayer 
to appear is the decrease of the density in the first layer by the promotion of the 
second layer. This mechanism was found in the computer simulation studies for 
methane films on graphite [63] and for argon films on (100) MgO [157]. Note 
that this last system is characterized by the square symmetry of the surface lat- 
tice. In the case of model systems discussed above, it was also found that sub- 
monolayer films melt without any effects due to the promotion of the second 
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layer, while the melting of a dense monolayer solid is accompanied by the 
transfer of particles from the first to the second layer (see Figure 16). 
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Figure 16. Temperature changes in the densities of the first and the second layers in the 
films of different total number densities (shown in the figure) formed on the (100) plane 
of fcc model crystal characterized by the corrugation parameter y, =0.8 and obtained 
from the canonical ensemble Monte Carlo simulation study. The upper (lower) family 
of curves corresponds to the first (second) layer. (From Ref. 68) 


Although the above discussion presented only some selected examples of 
experimental and theoretical studies of ordering phenomena in adsorbed films, 
it demonstrated that the structure and the thermodynamic properties of adsorp- 
tion systems are very sensitive to the details of the molecular interactions oper- 
ating in such systems. In particular, we tried to show that one of the most im- 
portant factors that determine the behavior of adsorbed layers is the corrugated 
surface potential and the symmetry of the surface lattice. Also, we attemped to 
show that computer simulation methods are powerful tools, which allow to get 
important new insights into the physics of adsorption systems as well as lead to 
the results allowing for a better understanding of several experimentally ob- 
served phenomena. 
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